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Abstract. We study inviscid limits of invariant measures for the 2D Stochastic Navier-Stokes equations. As 
shown in [Kuk()4] the noise scaling \pv is the only one which leads to non-trivial limiting measures, which are 
invariant for the 2D Euler equations. We show that any limiting measure /io is in fact supported on bounded 
vorticities. Relationships of to the long term dynamics of Euler in the L°° with the weak* topology are 
discussed. In view of the Batchelor-Krainchnan 2D turbulence theory, we also consider inviscid limits for the 
weakly damped stochastic Navier-Stokes equation. In this setting we show that only an order zero noise (i.e. 
the noise scaling v") leads to a nontrivial limiting measure in the inviscid limit. 



1. Introduction 

We consider incompressible Euler's equations and the randomly forced incompressible Navier-Stokes 
equation on a two-dimensional torus = IR^/A, where A c is a lattice. ' The Euler equation will mostly 
be considered in the vorticity form 

5fW-HM- Vaj = 0, (1.1) 

where we assume that w(x, t]dx-{i and the velocity field w = [u\, U2) is determined by o) from the 
equations 

curlM = V"*" • M = 3iU2 - 52^1 = ai, divM = 0, / u = 0. (1.2) 

In terms of the stream function, u = V-^y/, Ai// = o), with the usual notation V-^y/ = (-Vx2> Vxi)- The Navier- 
Stokes equation will be written either in the velocity formulation 

dtU + uVu-\ vl\u = f, (1.3) 

P 

where p is the (constant) density and the forcing term / = f{x, t) satisfies fj2 fix, t)dx = for each t, or in 
the vorticity formulation 

5t(U + M- Vw - vAw = curl/ , (1.4) 

with the relation between u and o) given as above. The particular stochastic form of / will be discussed 
below (1.10). 

In this paper we study certain classes of invariant measures for (1.1) and related equations. In particular, 
we address regularity properties and relations to the long term dynamics of 2D Euler for a particular class of 
invariant measures which arise as an inviscid limit of the stochastic Navier-Stokes equations, with a suitable 
scaling of the noise coefficients. 



Date: February 5, 2013 

^The reason why we consider general flat tori, rather then just K^/Z^ is that the geometry of the torus might have some influ- 
ence on various predictions concerning the long-time behavior of the solution. This issue will however come up only tangentially, 
and it will not be important for the results proved in the paper. The reader can take - U^/Z^ or - U^IIitZ^ most of the time. 
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To discuss topics wiiicii come up in various accounts of 2D turbulence, we will also use a linear damping 
operator Y defined in the Fourier coordinates by 

Yuk = rkUk, (1-6) 

where jk ^ 0- Often it is assumed that jk 7^ only for a few lowest modes, but other options, such as 
Yu = ju for some 7 > (corresponding to = 7 for each k) are also possible. We will denote the operator 
Y with 7fc 7^ only for a few low modes by Fiow - Its precise form will not be important for our discussion. 

1.1. Two dimensional turbulence. The standard theory of 2D turbulence conjecturally describes the 
behavior of solutions of 

(9fM + M-VMH vAm+FiowM = /i (1-7) 

p 

where / = /(x) is a "sufficiently generic" smooth vector field, which is on the Fourier side supported in 
a few relatively low Fourier modes. One can expect that for many / the system will become "turbulent" 
for sufficiently low v, although it is important to keep in mind that there are examples where this is not the 
case, see [Mar87, CR07]. In the turbulent regime one expects the famous downward cascade of energy 
together with an upward cascade of vorticity, as conjectured by Kraichnan [Kra67] and Batchelor [Bat69]. 
this is why we need the operator Flow ■ In dimension n = 3 the operator Flow should not be needed and we 
expect the so-called Kolmogorov-Richardson energy cascade. The striking feature of these phenomena is 
that, conjecturally, as v — ► 0+, the velocity field u should satisfy some bounds independent of v, such as 

II M II i<» < C, C independent of v . (1.8) 

Moreover, there are conjectures as to how energy will be distributed in the Fourier modes (see, for instance, 
the classical works [KM80, MWC92, Fri95, TsiOl, FJMR02, Tab02]). Rigorous treatment of these sce- 
narios seem to be our of reach of the present-day techniques and we have nothing new to say in this direction. 
Note however that the following ID model given by the Burgers equation 

d,u^uu.-vu.. = m, .eWZ, fuix,tUx=0, f fMd. =0, (1.9) 

is treatable. The behavior of the solutions for v — ► 0+ can be studied in detail via the Cole-Hopf transforma- 
tion. In particular, the bound (1.8) can be established in this case. 

Instead of relying on the chaos produced by the presumably complicated dynamics of (1.7) for low v, 
we can input "genericity" into the system by considering a "random" /. This point of view may be traced 
back to Novikov [Nov65] (see also, e.g. [BT73, VKF79]). We may take for example 

f[x,t) = aY,bkekix]W''it), (1.10) 

k 

where ej^ix) is given by (1.5), the sum is finite, over a few relatively low modes, a is a suitable constant 
of order 1. The are independent copies of the standard Wiener process (Brownian motion) so that 
are white noise processes and hence are stationary in time. After a suitable non-dimensionalization, a 
representative case of (1.10) is when l^fcP and a are both of order unity."^ 

^ The Fourier representation in set up in the following way. For scalar functions v on we will write vix] = Y. keZizA' vj^e'^'" , 
where A* is the lattice dual to A. Our functions v will satisfy fj2 vix) dx = which is the same as vq = and the above sum will 
always be taken over the non-zero elements of ZttA* . If there is no danger of confusion we will write simply vix) = Xfc vjce^^^. 
Divergence-free vector fields / will be written as fix) = Y.k fk^kM, where 

,,W = [zi^.!^l,«-*=^. (1.5) 



\k\ \k\ 

We note that in our normalization we have ^fr^ Jj2 I dx = Y.k I "fcl^- pf!]" /t^ I/I^ '^^ - ^k l/fcl^- 

The assumption of "sufficient generiticity" is important, as we need that the system creates enough "chaos" for low v. 
^If we wish to consider dimensions of the various quantities, one natural choice seems to be to take ej; dimension-less, bj^ of 
the same dimension as u (which is [length] [time]"^), and a of dimension [tinie]~2 , so that aW^it) has dimension [time]"'^. 
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With a random forcing of the form (1.10), the equation (1.7) can then be viewed as a stochastic equation. 
With some additional assumptions, there will exist a unique invariant measure fJ. = fJ-v for the process defined 
by (1.7), see e.g. [FM95, DPZ96, Mat99, Mat02, BKLOl, KSOl, KS02, Mat03, MP06, HM06, KuplO, 
HMll, Debll, KS12] and containing references. Relations conjectured by Kraichnan's theory would then 
be satisfied in a suitable mean value sense. The benefit of working with the random forcing is that even 
though we still cannot make much progress on establishing Kraichan's conjectures in this setup, we now 
at least have a quite canonical object for our analysis, the measure ^.^ Indeed, the above mentioned works 
establish ergodic or even mixing properties of /i. These properties provide some theoretical justification for 
the measurement of the physical quantities described in turbulence theories. In the deterministic case the 
measure // should presumably by replaced by a suitable invariant measure on the attractor, see e.g. [CF88, 
FMRTOl]. 

One can of course go through similar considerations in three dimensions, but in that case the lack of 
rigorous results concerning the basic existence and uniqueness questions for the Navier-Stokes solutions 
prevents obtaining rigorously the measure jj. above (or its analogues on the purported attractors in the deter- 
ministic case). Note however the recent works [DPD03, FR08, Debll] on weaker notions of solutions and 
associated invariant measures. 

1.2. Kuksin measures. Kuksin [Kuk04], see also [KP05, Kuk06b, Kuk06a, Kuk07, KukOS, Shill, 
KS12], suggested to study of a different limiting regime and put in (1.10) 

a = c\/v, (1.11) 

where c is a constant independent of v.^ Assume we have this scaling and omit the terms u ■ Vtt, Yjow from 
the equation. Then for each Fourier mode we have 



|fc|2v \k\^' 



H\Ukn = ^-^^ = \bk\'—,, (1.12) 



a bound independent of v. In the non-linear case a similar bound can be obtained, see [Kuk04, KS12] 
and (1.17). With scaling (1.11) the operator Y is no longer needed^ and the invariant measures ii = fiy of the 
process generated by the equation 



du+ I M- VuH vAm 

p 



dt=cVvY,bkek[x]dw''it] (1.13) 

k 



have a meaningful limit (perhaps after passing to a suitable subsequence). At the level of the vorticity 
0) = curl u we have 

dcD + {u-ya)-v\a})dt = c\/vY,gkdw''it)e''''', gk = \k\bk, (1.14) 

fc 

which is the form which we will mostly work with.^ A deterministic version of the situation considered by 
Kuksin would correspond to setting 

f=cvf (1.15) 

in (1.3). This situation is relevant for Section 5. 

Let jUy be an invariant measure on the space of the vorticity functions o) of the process defined by (1 . 14). 
Note that for a sufficiently fast decay in the fo^'s in (1.13), these measures are supported on smooth 



^The case of the Burgers equation (1.9) with stochastic forcing can be analyzed rigorously, see for example [EKMSOO]. 
^If our quantities are not dimensionless and we use the same dimension count as in a previous footnote, then c should have 
dimension [length]"^. 

-7 

Indeed in this scaling, (1.11), the term Y leads to a trivial limiting measure ^ as we establish rigorously in Section 6 below. 
^The term -vAw can be replaced by more general dissipation, such as fractional Laplacian v(-A)"; see Remark 4.5 below. 
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functions, see e.g. [KS12]. Also, under rather general assumptions on these bj^'s such a measure is unique, 
cf. references above. Applying Ito's formula 

d j -0)^ dx = j [a)da) + ^dtodo)] dx (1-16) 

and taking the expectation, we obtain 

E[j^JVoj|2rfxj= j \\Va)\\l,dfiy{a)) = ^Y.\Sk\^- (l-l^) 

Due to this bound, as v — ► 0+ , the family of measures /iy has a subsequence converging weakly^ to a limit 
/Uo, which is a measure supported in the space 

= i^a)eH\j^), j ^cjdx = 0^, (1.18) 

with ^ 

E(||Va;||2j= [ \\ya)\\l,dno{a)]<^Y.\Skf- (1-19) 
J ^ k 

Kuksin [Kuk04] proves that this measure invariant for the evolution given by Euler's equation (1.1). See 
also Kuksin and Shirikian [KS12] and references therein for many other interesting properties, such as the 
non-triviality of the measure. We will call the measures constructed in this way Kuksin measures. 

1.3. Main Result on Kuksin measures and the dynamical systems approach to 2D Euler. We now 

show that the Kuksin measures are closely related to the long-time behavior of solutions of Euler's equation. 
One of our main results in this paper will be the following: 

Theorem 1.1 (Kuksin measures are supported on L°°). Let jibe a Kuksin measure as above. Then 

\(i)\\ioo d^Qicj) <+oo. (1.20) 



We shall discuss the outline of the proof of Theorem 1 . 1 in Subsection 1 .4 below. A detailed statement of 
the result and its proof is found in Section 4. 

In particular, from (1 .20) we see that ^ is supported on L°°. This is important, as the space L°° is prob- 
ably the most natural space (for the vorticities) in which to consider the 2D Euler equation when studying 
the long-time behavior of the solution. This is due to the following facts: 

(i) The initial value problem for (1 . 1) is well-posed for in a classical result by Yudovich [Jud63]. 

(ii) Let R>Q and let X = Xr = [cje L°°, fj2 w = 0, lloiHioo < R}. Equipped with the weak* topology, 
the set X is a compact metric space, which we will denote by (X, w*). One can check that the 
proof of Yudovich's theorem actually give a stronger result: namely, the Euler equation (1.1) gives 
a well-defined dynamical system on X (for any R> 0). A proof of Yudovich's theorem which 
can be easily adapted to prove our statement here can be found for example in [MB02]. 

From Theorem 1 . 1 we hence see that Kuksin measures (restricted to X) give natural invariant measures for 
the Euler evolution on X. The functions on which the measures are supported have additional //q -regularity. 

Note that one can also construct non-trivial measures on X which are invariant under the Euler equation 
directly: we know that the energy functional'*' 



1 r 1 1 r 1 

= — — I -\u\^dx= — — / — wcjdx (1-21) 

is continuous on (X, w*). Therefore the energy level sets Xe = Xr^e given by {to e X,S(_a)) = E} are compact 
subsets in X which are invariant under the Euler equation (due to the energy conservation). By the classical 



More precisely we have that /^a f[(t))djivj (w) converges to Jji f{<i))dfiQ[(i)) for every continuous, bounded real valued test 
function /. In fact the convergence holds also in H^~^ for any e positive. 
^'^More precisely, energy per unit mass. 
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Kryloff-Bogoliouboff procedure, every non-empty Xr^e supports an invariant measure. This measure cannot 
be trivial when E>0 i.e. supported at to = 0. There are additional conserved quantities for the evolution by 
Euler's equation, namely the integrals 



but these quantities are not continuous on [X, w*), making their implications for the dynamics on [X, w*) 
more subtle. 

In Section 2 we discuss various hypothesis for the long term dynamics of the (1.1) on the phase space 
{X, w*). These hypothesis may illuminate further possible structure of the support of the invariant measures 
//o- Two extreme scenarios present themselves. On, the one hand, taking the view of statistical mechanics, 
we may predict long time behavior from the maximization of various notions of "entropy" subject to the 
constrains of the Eulerian dynamics. In many cases these "entropy maximizers" may have a fairly simple 
shear flow like structure. Thus in this scenario we would expect that the Kuksin measures would be sup- 
ported on steady flows with a relatively simple topology. At the other extreme we might suppose that all 
of the solution trajectories of the Euler dynamical system are pre-compact in L^. In this case many of the 
Casimirs (1 .22) must be conserved at the end-states which would suggest that /io has a much richer structure. 

There is some evidence for both of the above scenarios. On the one hand we do not have a single 
example where it is proved that an initial condition yields an orbit which is not precompact in L^. Moreover 
a recent result in [Svel2] (which we recall in Theorem 2.1 below) rigorously shows that at least some 
such precompact orbits must exist. On the other hand recent numerical result of [BS09] suggest that is 
concentrated on certain laminar states obtained as an "Entropy maximization". It seems unlikely that either 
of these scenarios holds universally and that the structure of /io is given by an intermediate situation. 

1.4. Moser iteration for SPDE and applications to L°° estimates for stationary solutions. We now 

turn to discuss some aspects of the proof of Theorem 1.1. We will see that the main ingredients involve 
a suitable rescaling of the equations and then developing a Moser iteration scheme for SPDEs of drift 
diffusion type which evidences a parabolic regularization from to L°°. The detailed proofs are given 
below in Section 4. 

A natural rescaling of time makes the interpretation of the measures jdy and the Kuksin measures //q 
perhaps more transparent. If we replace the function u(x, t) by u[x, t) = u[x, tlv) and replace W^[t) by the 
equivalent process W(_t) = s/vW^it/v), we obtain, after dropping the tildes we obtain 



Note that the measure Hv is also an invariant measure for this process. See Section 2.1 below for some 
motivating discussion of analogous finite dimensional situations. 

As V — in ( 1 .23) the drift velocity v~^u grows unboundedly. As such to obtain Theorem 1 . 1 need find 
a way estimate L°° norms of solutions to equations of the form 



with constants that do not depend on the size of the sufficiently smooth divergence free drift velocity a. 

In the deterministic case, one usually obtains such drift-independent L°° bounds either by appealing 
to maximum principle-type arguments, or by using estimates, with p independent bounds, and passing 
p — ' oo. Neither of these direct approaches appear to be available in the stochastic case. The first approach 
seems to fail since one cannot exchange E and sup^, and due to the lack of smoothness in time of the 



Since the noise in (1.24) is additive, one could shift the equation by subtracting the solution of an Omstein-Uhlenbeck 
process, but the L°° bounds one obtains on the resulting random PDE appear to depend essentially on the size of the drift velocity. 




(1.22) 




(1.23) 



dct)+ [a- Voj- Aa))dT = cj^gjcdw'^, V-fl = 

k 



(1.24) 
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Stochastic terms in (1.24). On the other hand, for Lf bounds, a direct application of the Ito lemma to (1 .24) 
yields 



p{Aa},oj\cij\f 



(p-2)/2|| 



\dt + pY,{g,„(o\a)\P-^)dw'' 



(1.25) 



where we have used that a is divergence-free. The Ito correction term in (1.25) grows quadratically in p, 
which is too fast. On the other hand, letting X = \\ci)\\^p, one may apply Ito's lemma to cpiX) = (1 + X)^'P/p 



(see e.g. [KrylO, Remark 5.2]) and prove using standard estimates that 

\\(0{t) 



sup 

p>2 



E sup ■ 

mo.T] 



||2 
II LP 



P 



<c 


( 

1 + sup 


/ 

E- 


1/, l|2 ^ 


+ vr||cr||Lcx= 




p>2 




P J 





where u) is the solution of (3.2). This however does not yield a bound on E|supfe[Q supp>2 ||a»(f)|||p/pj. 

Since (1.24) is a parabolic SPDE, in the spirit the classical DeGiorgi-Nash-Moser [DG57, Nas58, 
Mos60] theory for deterministic parabolic PDEs, one may expect an instant regularization of the solution. 
The difficulty in carrying over this program lies in treating the stochastic forcing term in ( 1 .24) and obtaining 
bounds which are independent of the size of drift velocity a. In the deterministic case, for drift velocities 
that are divergence free, one obtains the to regularization of solutions to the parabolic equation, with 
bounds that are independent of the drift using e.g. the elegant argument of Nash [Nas58]. Drift indepen- 
dent bounds for a deterministic analogue of (1.24) have also been obtained using Moser iteration, see, e.g. 
[Kuk99]. Therefore, one may expect that the same result holds for stochastic drift-diffusion equations such 
as (1.24). 

In order to treat the stochastic term, it turns out that the iteration technique introduced by Moser [Mos60] 
is better suited in view of the Ito formula (1.25). This fact was recently observed in the context of 
semilinear SPDE in [DMS05, DMS09] where the authors obtain an LP° maximum principle. We however 
cannot appeal to these results since they rely essentially on the fact that the initial data already lies in L°°. 
By ( 1 . 17) we only have v-independent bounds on the statistically steady solutions of ( 1 . 14). To overcome 
this difficulty we prove in Theorem 4.1 (see also Remark 4.4 below) that the solution oj{t) of (1.24) lies in 
L°° (in x) for arbitrarily small positive time t: 



E sup 

tElT,2T] 



\a){t)\\Loo<C[l + T-^'*)E 



l^^llLi([0,2ri;L2) VC^I 
k 



(1.26) 



for r < 1/8, where the constant C is independent on a. To the best of our knowledge the parabolic regular- 
ization estimate (1 .26) is new in the context of SPDE. As in [DMS05, DMS09], one of the main differences 
between the stochastic Moser iteration (see the proof of Theorem 4. 1) and the classical approach for deter- 
ministic PDE naturally arises due to the random forcing. In view of the Burkholder-Davis-Gundy inequality 
we need to bound quadratic variations of the Martingale on the right side of (1.25), and hence the integra- 
bility in time needs to be twice that in space in order to close the iteration scheme (cf. (4.32) below). 

In view of the predictions made by Statistical Mechanics arguments regarding the "end states" of the 2D 
Euler dynamics, and having already established that the Kuksin measures are supported on (T^) nL°°(T^), 
we beheve that: 

Conjecture 1.2. Kuksin measures are in fact supported on continuous vorticities. 

The immediate difficulty which arises in proving this conjecture is that even in the deterministic case, 
for the two-dimensional linear parabolic equation 



dtV + bix,t)-Vv-\v = f, y-b = 0, 



(1.27) 



the size of the smooth drift comes into play for the DeGiorgi-Nash-Moser proof of Holder regularity. If the 
drift is rough, one may even construct solutions that are not continuous functions for all time, although they 
obey the L°° maximum principle (see, e.g. [SVZ12]). 
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On the other hand, one of the key ingredients of the proof of Theorem 4. 1 was the (statistical) stationarity 
of the solution to (1 .23). As a deterministic toy problem one may hence consider time-independent solutions 
of (1 .27), with drift b{x). In this case, following the ideas in [SSSZ12] and an elliptic Moser iteration we are 
indeed able to prove in Theorem 5.1 below that the solution obeys a drift-independent logarithmic modulus 
of continuity, and is hence uniformly continuous. The analogy between time-independent solutions to (1 .27) 
and statistically stationary solutions of (1.24) is however tentative at best. 

1.5. Inviscid limits for damped models; different scalings. In view of the foregoing discussion con- 
cerning the Batchelor-Kraichnan theory of 2D turbulence, it appears that when working on the periodic box 
the stochastic Navier-Stokes equations should be augmented (as in (1.7)) with a suitable damping term Y 
to prevent a pile up of energy at large scales. Note that the Y term also frequently appears in geophysical 
models closely related to the 2D Navier-Stokes equations to account for friction with boundaries. In these 
situations, if the scaling in this damping term is held fixed as v ^ 0, then a different scaling must be in- 
troduced for the noise in order to obtain a non-trivial inviscid limit in the class of the associated invariant 
measures. 

To this end, we consider operators of the form Y = Yj^ = tA"^ = t(-A)"T''^ and t > 0, 7 e [0, 1) and 
study weakly damped and driven stochastic Navier-Stokes equations of the form 

du + {Yu+u-Vu + Vn-vAu)dt = v"pdW, V-m = 0, (1.28) 

for different values of a e D?.'" As above for the undamped case (1.17) energetic considerations allow us to 
deduce the correct scaling with a in (1.28). Consider a collection of invariant measures {;u"}v>o for (1.28). 
Let be stationary solutions of (1.28) corresponding to /z" and denote = V-*- • m^. Applying the Ito 
lemma to the vorticity formulation of (1.28) and using stationarity one deduces that: 

E(v||Va;^|li. + ||Fi'V|li.)~ll-llL- 

Making use of the above relation, we will show below in Theorem 6. 1 that a = is the only relevant scaling 
for (1.28). Here stationary solutions of a damped stochastic Euler equation arise. See also [BF12, Bes08]. 

Organization of the Paper. In Section 2 we review some notions related to the time- asymptotic behav- 
ior of the 2D Euler equations. Our discussions in this section allow us to make some hypotheses regarding 
the structure of the support of Kuksin measures in this context. Section 3 recalls the mathematical framework 
for the Navier-Stokes Equations and its associated Markovian semigroup. We also review some properties 
of Kuksin measures established in previous works. Section 4 is devoted to the proof of the main theorem. 
Here we detail the Moser iteration scheme which addresses a more general class of drift-diffusion equa- 
tions. Section 5 concerns a deterministic toy model for stationary solutions of the stochastic Navier-Stokes 
equations. We establish a modulus of continuity for this system. The final Section 6 we consider a weakly 
damped stochastic Navier-Stokes equation and establish inviscid limits in the appropriate scaling for this 
model. 

2. Long term behavior of 2D Euler and related systems; connections to invariant measures 

In this section we discuss some aspects of the long time dynamics in (L°°, w*) of solutions to 2D Euler 
and the relation of this behavior to possible properties of the Kuksin measures, which are now accessible 
due to Theorem 1.3. We begin with some motivation from finite dimensional Hamiltonian systems. 

2.1. Noise scaling limits in finite dimensions. A finite dimensional situation related to the above is 
studied in the theory of the small random perturbations of dynamical systems. Let 

x=b{x) (2.1) 



As explained above it is of interest to consider Y acting only at the largest scales; cf. (1.6). This situation is more delicate 
to analyze rigorously and would seem to require the establishment of suitable "hypocoercivity" properties for (1.7). 
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be a dynamical system in IR". Consider its stochastic perturbation 

dx=b{x)dt+\/eQdW, (2.2) 

where W = {W^,...,W"-) are normalized independent Wiener processes and Q is a matrix. By setting 
xit) = xiet), W(_t) = -^W{£t] and t = et and dropping the tildes, we obtain 

dx=-b{x)dT + QdW. (2.3) 

£ 

Such systems have been extensively studied, see e.g. [FW12]. In case of measure -preserving flows, the 
behavior of (2.3) as e ^ 0+ can be understood from the following picture: under some assumptions the 
equation x = ^b[x) takes the trajectories very quickly through "ergodic components", and hence for £ — ' 0+ 
the system (2.3) should in some sense describe a diffusion in the space of the ergodic components. 

Equation (1.13) (or its rescaled version (1.23)) is of a slightly different nature that the perturbation of 
Hamiltonian systems considered in [FW12], in that we add not only a small noise, but also small damping. 
Such procedure can be illustrated by a simple example of the Langevin oscillator: 

Example (The Langevin oscillator). We consider a simple Id harmonic oscillator with damping and 
random forcing 

mq + Yq + Kq = aw , (2.4) 

where is a normalized Wiener process and m,Y,K> 0. Letting p = mq as usual, it is easy to check that 
the (unique) invariant measure of the system 

^. r ■ (2-5) 

p = -xq- j^p + aw 



is given by the Gibbs measure 
where the Hamiltonian H given by 



dp=-e-^^^l-P^dqdp, (2.6) 



2m 2 

and Z = Z{fi,K,m) is a suitable normalizing constant. We see that from the point of view of Statistical 
Mechanics the quantity y corresponds to (a multiple of) temperature. A similar calculation can be done for 
a general one dimentional Hamiltonian of the form 

H=^ + V{q). (2.8) 
2m 

In higher dimensions one can also calculate further examples; especially when the Hamiltonian is qua- 
dratic. The invariant measure does not necessarily have to be the Gibbs canonical measure as in (2.6). If the 
damping and the forcing are taken to with the analogue of the ratio y converging to a limit, the invariant 
measure will converge to an invariant measure of the original Hamiltonian system. For example, in the case 
of a completely integrable n-dimensional system with the full system of mutually commuting integrals of 
motion /i , . . . , /„ the limiting invariant measure can be expected to be of the form 

dp = ^^'^'■^^ ^"^ dqi...dqndpi,...,dp„, (2.9) 

where the function (p will depend on specific choices of damping/forcing. We see that the vanishing damp- 
ing/random forcing method can be viewed as a way of producing suitable statistical "ensembles", closely 
related to those used in Statistical Mechanics. Considerations in this direction in the context of the KdV 
equation can be found in [Kuk07]. In the terms of Statistical Mechanics the ensembles produced by this 
method are related to "canonical ensembles". One can also consider the "micro-canonical ensembles", 
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which in the last examples would simply be given (under some "genericity assumptions") by invariant mea- 
sures on the tori 

/l = Ci, /2 = C2, ...,/„ = c„. (2.10) 
Under suitable assumptions, invariant measures on these tori are the "irreducible components" of the mea- 
sures (2.9). By analogy, we see that Kuksin measures should be related to the Statistical Mechanics of 
Euler's equation. Their decomposition into irreducible components should give "ergodic components" of 
the Euler evolution. However, this analogy may break down due to infinite dimensional effects. As v — ► 0+, 
it is conceivable that fast Euler evolution moves enstrophy to high spatial frequencies (in the Fourier spec- 
trum), so that "complexity" is lost (by disappearing to infinity in the Fourier space) and Kuksin measures 
may conceivably be supported on some relatively simple sets, perhaps even equilibria. This would be an 
infinite-dimensional effect,' ' which does not have an analogy for finite dimensional or completely integrable 
hamiltonian systems. This is discussed in more details below, but still at a heuristic level. We were not able 
to obtain rigorous results in this direction. 

2.2. Long-time behavior of solutions of Euler's equation. Equation (1.23) together with some anal- 
ysis of the long-time behavior of solutions of Euler's equation seems to give some good hints about what 
one should expect concerning some of the properties of Kuksin measures. We recall some of the expected 
properties of the Euler solutions. 

We consider equation (1.1), with the conventions (1.2). We also recall the obvious identity 

/ LDix,t)dx = 0, teU. (2.11) 

Let us start with some classical observations about the long-time behavior of the solutions of (1. 1) starting 
from initial data oiq e L°°, with ||wolli°° = R- Let X = Xr be the ball of radius R in L°°, taken with the 
weak* topology. In addition, we can impose the constraint fj2 o) dx = on the functions in X. The space 
(X, w*) is a metric space and, as already discussed, the Euler evolution (1.1) gives a well-defined dynamical 
system on X. We denote the H-limit sets by 



Q+ =n+(wo) = nt>o{w(5), 5> r,w(0) = wo} . (2.12) 

We also introduce 

^ct)o = {^^0 °h, h : ^ is a volume-preserving C^-homeomorphism} (2. 13) 

and 

0<^o,E = (^a>on{(^,<Si(^) = E}. (2.14) 

It is not hard to see that the weak* closure of 0,^^ , denoted by is a closed convex subset of for any 
p > 1. Letting S{a)o) = E, we clearly have 

n^^^X.E- (2-15) 

There are various conjectures concerning the long-time behavior of Euler solutions motivated by the notion 
of "mixing", see [Mil90, Rob91, Shn93, Svel2]. We can think of the fluid as consisting of fluid particles, 
with each fluid particle having a fixed value of vorticity permanently attached to it. The fluid motion mixes 
these particles, with the vorticity remaining attached to each particle. The most naive conjecture could be 
that for large times the vorticity is everywhere mixed, corresponding to the weak* convergence of w{f) to 
as t—'oo. This would mean'"^ 

Ci+ = {0}. (2.16) 
In the Fourier space this would correspond to the movement of all energy''' of the solution towards 
larger and larger frequencies as f — ► oo. If £" 0, then (2. 15) provides an obstacle to this. The energy cannot 



1 ^ 

The effect is closely related to Landau damping, see for example [MVll]. 
^'^This presumably happens if we consider the Burgers equation with the scaling (1.13). It should not be hard to verify that for 
the Burgers equation the Kuksin measures are trivial. 

'^It would be more precise to say enstrophy, but in the situation here this does not make a difference. 
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all move to high (spatial) frequencies, as the energy functional S is weakly* continuous. We can "fix" this 
by trying to "move" as much as possible energy to high frequencies which is still compatible with (2.15). 
More specifically, we can try to minimize 



|T2|Jt2 Y 

subject to the constraint 

^^e^r„'£- (2-18) 
Note that /(w) is preserved during the evolution, but it can conceivably drop on the "end-states" Q+, as it is 
not weakly* continuous. 

Minimizing / (subject to (2.18)) is of course the same as maximizing -/(oi) subject to (2.18). More 
generally, we can take a concave function F and maximize /f (w) given by (1.22), subject to (2.18). The 
quantify /f (w) can be called the entropy of the "configuration" di and the above principle is then nothing 
but the usual entropy maximization under given constraints, as well-known from Statistical Mechanics. The 
entropy /f is could be considered as too simple, the usual entropy in Statistical Mechanics is based on 
suitable "counting of states". Closely related to the notion of entropy is A. Shnirelman's notion of mixing 
in [Shn93]. 

There are indeed more sophisticated notions of entropy, see for example [Mil90, Rob91, Tur99, Svel2], 
which can be more "non-local" than the above. For example, let wq = aiXA], where A; is a division 
of into disjoint measurable sets with | A/| = x/IT^I, and L/ ^J/x/ = 0. Then, the closure of defined in 
(2.13) is 



-w* 

Wo 



and one can define the entropy (generated by wq) as 

S(w) = S(^„(w) = sup i f ^-p/logp/rfx: a»(x) = ^a/p/(x), 0<p/<l, Y.P^^^\- (2.20) 

\\^ \J^'^ I I I \ 

This entropy leads to the theories of Miller and Robert, [Mil90, Rob91]. When the division = U/A/ has 
only two sets A\ and A2, then this entropy is of the form If for a suitable F. For example, when we only 
have two sets and a\ = -a2-\, then 



F(W): 



1 +ai 



1 + W 



log ^|- ^ log i^l. (2.21) 



The maximizers of the entropy subject to given constraints are steady-state solutions of Euler's equations 
given by stream functions satisfying 

= H{y/) (2.22) 
for a suitable function H depending on Ip. These should be the "end-states" of the evolution if the actual 
evolution and Statistical Mechanics lead to the same conclusions. 

As is usually the case with predictions based on Statistical Mechanics considerations, it is difficult to 
decide whether the actual dynamics of the equation produces the behavior expected from entropy maxi- 
mization, assuming all known the constraints have been taken into account. In fact, we do not have a single 
example which in which it would be rigorously established that the trajectory 

Qf = u{w(5), s> /■} (2.23) 

is not pre-compact in (and hence any for p e [l,oo)). On the other hand, it is useful to recall the 
following result from [Svel2] : 

Theorem 2. 1 (Existence of precompact orbits [Svel2]). The omega-limit set n+ of any trajectory 
always contains an element (Oq whose trajectory is pre-compact in L^. 



^''The proof of this statement is very simple: maximize some entropy If with a strictly concave F over Q.+ . 
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2.3. Possible Consequences for Kuksin Measures. In view of the above discussions concerning the 
long term dynamics of 2D Euler we now introduce two extreme scenarios: 

Scenario A: Euler solutions weakly* approach entropy maximizers. Let us assume that our torus is 
= IR^ / flZ ® foZ with 0< a<b. Let us further assume that 

(i) All entropy maximizers for Euler solutions (with given constraints) are shear flows independent 
of Xi . This is in fact not very far-fetched. It has been established rigorously for sufficiently small 
energies and local entropies Ip with F strictly concave. See [FS]. 

(ii) All solutions weakly* approach these shear flows as ? — ► oo. This would of course be a very strong 
statement which we do not really expect to be true. However, if (i) is correct, then the Statistical 
Mechanics predictions would suggest exactly this conclusion. 

If this scenario holds, then one can expect the Kuksin measures to be supported on the steady-state shear 
flows. Indeed, from the re-scaled equation (1.23) we see that as v ^ 0, the fast Euler dynamics will drive 
the solution towards the shear flows, whereas the term Aai will be quickly damping the high frequency 
components of o) generated by the Euler evolution. This scenario is genuinely infinite-dimensional: all the 
complexity of the Euler dynamics and the initial data will disappear into the high frequencies, and will never 
"return". Such behavior does not have an analogy in finite-dimensional systems or in completely integrable 
systems. 

Scenario B: all solution orbits are pre-compact. Let us assume that the solution trajectories Qf 
in (2.23) are pre-compact in L^. This may be unlikely, but it has not been ruled out. In this case the 
weak closures of these trajectories will be the same as the strong closures and all the functionals Ip will 
be conserved on the "end-states". In particular, the mixing envisaged by the statistical mechanics approach 
will never take place. In this case the Kuksin measures will have much richer structure. Their "irreducible 
components", similar to the measures on the tori (2.10) in the example leading to (2.9), will be supported 
on the closures of the -compact trajectories, which will play a role somewhat similar to the ergodic 
components in finite-dimensional Hamiltonian systems. In this scenario many features familiar from finite 
dimensions or completely integrable systems will still be present. 

We conjecture that neither of these scenarios is quite true, but that the real behavior will be intermediate 
between these two extremes: non-trivial L^-precompact trajectories will exist, but initial data leading to 
them will not be "generic". The Kuksin measures will be supported on such trajectories. Depending on 
our degree of optimism, we can hope that these solutions represent a type of a weak* attractor for all Euler 
solutions. 

3. Some results regarding invariant measures and inviscid limits 

In this section we first recall some elements of the mathematical analysis of the stochastic Navier-Stokes 
equations and its associated ergodic properties. This allows us then to summarize some of the analyti- 
cal properties enjoyed by the Kuksin measures, established in previous works (cf. [KS12] and references 
therein). 

3.1. Mathematical setting: stochastic 2D NSE and its Markov semigroup. We consider the Sto- 
chastic Navier-Stokes Equations on a periodic box 

du+[u-Vu + y7i-v\u)dt=\/vpdW=\/vY,Pkdw'^, V-u = 0, u[0] = uq. (3.1) 

k 

As discussed in the Introduction, in order to consider the inviscid limit v — ► 0, we use the scaling \/v for the 
noise coefficient. Typically we will use the vorticity formulation of (3.1). Taking o) = V"*- • m we obtain 

da) + iu-ya)-v\a))dt=\/vadW=\/vY,ajcdw'', a)(0) = too. (3.2) 

k 

We will assume that fj2 toodx = and fj2 crdx = 0, which implies that the solution is always mean zero. 
Note that u can be recovered from oj via the Biot-Savart law. 
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Let us now set some notation used throughout the rest of the work. We denote the Sobolev spaces 



with = Lpg^ and take the usual norms and inner products donated by || ■ || (•, Ofc- We will denote the L^, 
p> 1 by II • IIlp. 

To emphasize dependence on initial conditions we will write O)^it,a)o) = w(t,a»o) for the solution of 
(3.2) with initial condition oiq- Assuming 

Y^\\ai\\l,<oo, (3.3) 

we have w(-,(yo) e C([0,oo);//°) nL^^^([0,oo);i?^) for any oje H'^. For k>l, assuming that 

Y,\\ai\\l,<oo, OA) 
I 

we also have the higher regularity properties for (3.2). If fc»o ^ H'^ then, for any to > 0, w(-,(0o) £ C([fo>oo);i?*^)n 
L2^J[fo,oo);//*^+i). Similarly, if e H'', (x){-,ojq) e C([0,oo);//'')nl2^^([0,oo);//*^+M. Note that the general 
well-posedness theory for the stochastic Navier-Stokes equations has been extensively developed. See e.g. 
[BT72, BT73, Vio76, Cru89, CG94, FG95, MR05, BreOO, BFOO, BPOO, MR04, GHZ09, DGHTll]. 
Notational Conventions for the Stochastic Terms: For brevity we will often write e.g. 

11^11^2 :=(X:ilo-/||2,)''' (3.5) 

when no confusion will arise from this abuse of notation. Similarly, for 2 < p < oo will also take 

i\ i/p 

j^^[Y,\aiix)fY'^dx (3.6) 

and, for p = oo, 

\\a\\L'^:=sup{Y^\aiMff'\ (3.7) 

We next recall some aspects of Markovian framework for (3.2). Take SSiH'^) to be the Borealian subsets 
of h'^. We define the transition functions 

Pf(wo,r) = P(w(f,wo)er) (3.8) 

for any f > 0, a»o e // and Y e SS{H^). Let CbiH'^) and Mi,{H^) be the set of all real valued bounded con- 
tinuous respectively Borel measurable functions on H^. For r > 0, define the Markov semigroup according 
to 

Pf0(wo) = E0(i^(?,wo))= f (l){(jj)Pt{o)a,da)) (3.9) 

which maps Mh[H^) into itself. Since cj[t,(t)o) depends continuously on a»o e H'^, it follows that Pt is Feller 
i.e. P[ maps CbiH'') into itself. Let Pr{H^) be the set of Borealian probability measures on h'^. Recall that 
jj. e Pr{H^) is an invariant measure for (3.2) if 

/ (p{(i)Q)dp[(i)Q)= \ Pt4){coo)diJ.{(Oo), for every f>0. (3.10) 

For further generalities of the ergodic theory of stochastic partial differential equations see e.g. [DPZ96, 
KS12]. 
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3.2. Existence and Uniqueness of invariant measures for SNSE. For each v > 0, there exists an 
invariant measure /iy in for (3.2). This can be established using the classical Kryloff-Bogoliouboff 
procedure, [KB37] by proving the tightness in H° of a sequence of time average measures starting from any 
convenient initial condition. Note that if (3.4) holds then it is not hard to show that /Zy is supported on H'^, 
[KS12]. We will denote by <w^{-) a statistically stationary solution of (3.2) associated to ^Uy. In other words 
/Uv(-) = P(w^(f)e-). 

The invariant measures /Zy along with the associated stationary solutions to^ satisfy the balance relation: 

IlmI. (3.11) 

^ k 

lL\\^k\\l2. (3.12) 

k 

We derive (3.11), (3.12) by applying the Ito formula to, respectively to (3.1), (3.2) for (o^. For example, 

E||(y(r,too)ll^2 + vEJ \\a}[s,a}o)f^ids = E\\a)o\\^^2 + vT\\a\\^^2 (3.13) 

One can also use the Ito formula to prove that 

Eexp(<5||wp|^2)<C<oo (3.14) 

for some 5 > and a constant C that is independent of v. See e.g. [KS12] and containing references. Note 
that 

ENsIIh^+i= / Noll^.«rfMv(t^o)<C(v) (3.15) 

where Civ) is finite. However, it is doubtful that we can bound this quantity C(v) independently of v for 
k>l. 

Remark 3.1 (Uniqueness of fiy for v > 0). For each v > the uniqueness of /iy is a much deeper 
question and requires the imposition of much specific conditions on a. One needs to establish smoothing 
properties of the Markov semigroup Pt (ellipticity or hypoellipticity of the Kolmogorov equation) and that 
a common state can be reach by the dynamics regardless of initial conditions (irreducibility). See, e.g. 
[FM95, DPZ96, Mat99, Mat02, BKLOl, KSOl, KS02, Mat03, MP06, HM06, KuplO, HMll, Debll, 
KS12]. Since the results we develop here related to inviscid limits do not require /iy to be unique, we do not 
impose such additional conditions on a. 

Remark 3.2 (Some explicit stationary solutions). We can identify some very special choices for a 
which allow us to obtain explicit statistically stationary solutions of (3.2). Suppose we have found any 
a)£ : ^ IR satisfying 

ue-Vme = 0, -Aai£ = Aoi^. (3.16) 

Here A > and ue is obtained from a)E via the Biot-Savart law. For example "laminar states" satisfy (3.16). 
Consider the process 

a)lit,x) = a)Eix)\/v f exp{-vX{t-s))dw}. (3.17) 

J-oo 

Let us note that to^ := ai£(x)zs(f), where zs is the unique stationary solution of the Id Omstein-Uhlenbeck 
(Langevin) process dz + vXzdt = \/vdW^. Here, zg is normally distributed with mean zero and variance 
(2A)"\ for each v > 0. Then is a stationary solution of 

da)+ iu-Vcj-vAa))dt= \/va)EdW^. (3.18) 

This may be checked, for example by using the mild formulation of (3.18). Hence, in this setting the 
invariant measure obtained as v ^ is also normally distributed around oje- 



E\\ci)l\\l2 = j WcDoWlidfivM = 



and 



E||w^"2 
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3.3. Previously established properties of Kuksin measures. The balance relation (3.12) implies that 
any collection of invariant measures = {|^v}v>o is tight and therefore weakly compact. We denote by 
jUo a limiting point of J?, and refer to these measures as Kuksin measures. Let us now recall some known 
properties of the measures ^o- We refer to [KS12] for the proofs of all the facts described in this subsection. 

Define = W^'^ (IR; //°) n i^^^dre, H^) and let J^e be the set of w e that satisfy the Euler equation 

dt(o + u-Va) = (3.19) 

in its vorticity form weakly for all t eU. Moreover if 0)1,0*2 £ J^e and wi(f) = 0*2 (?) for some t eU then 
(1)1 = 0)2- This follows from the methods of Yudovich (see [MB02]). Define n:^^Hhy nio)) = a){0), the 
fiber at f = 0. Note that, from uniqueness, it follows that n is injective on and let 

3f: = n[JfE)- 

It then holds that HoiS^) = 1. 

The Euler equation is well defined as a dynamical system on Indeed, for e there exists a 
unique solution a)[-,a)o) e c C(IR, //°) n L^^^{IR, H^). Indeed for Z: e IR define St : ^ , via Sf^o = 
(i)[t,a)o). We endow with the topology inherited from Jif, i.e. we take 

l + (g^Ar(w^,a)^) 

where 

Ilw(-,wj)-w(-,w^)ll2^^f. (3.21) 

We then have that {Sr}f£K is a group of homeomorphisms on Moreover ^0 is invariant for {St}t>o, i.e. 
^o(£) = ^oiS-\t)E] = iJ,oiS{-t]E], for all f e R and any E e SSiS^). 

Using local time techniques it may be shown that ^0 is "non-trivial" in the sense that it contains no 
atoms. In other words, for any o) e SC , we have that = 0. Further properties such as spacial ho- 

mogeneity, and higher moment bounds, and pointwise in space moment bounds for the measures /io are 
discussed in [KS12]. 

Remark 3.3 (Vortex patches and jUo). Let us observe that there are no vortex patch solutions in the 
support of /iQ. Indeed, for any open set O c IR^ let xo be the indictor function on O. Define 

= {xo '■ O cc T^, bounded, simply connected with smooth boundary 30} . 

For any ojo e 9, according to e.g. [MB02] there exists ?o > such that a)(f,a»o) e S^, for all f e [- fo. fol- As 
such, since //^ nS* = 0, for any wo e S^, w°(-,wo) t L^^^dR,//^) and hence w"(-,wo) t J^e- Thus Hai^) = 
since 

g^r)3^ = 0. (3.22) 



4. Invariant Measure Supported on L°° and Related Estimates 

In this section we establish uniform in v bounds on E||w^||i«>, where cjy are stationary solutions of (3.2). 
Our approach makes use of the Moser iteration technique and draws on earlier works in this direction in 
[DMS05, DMS09]. However we obtain parabolic regularization and time decay for the initial data compo- 
nent, which was not addressed in the above works. Let us note that in the deterministic case L°° bounds may 
be obtained without appealing to the Moser iteration: one carries out LP estimates, which take advantage of 
cancellation in the nonlinearity, and then sends p — 00. In the stochastic case, the Ito correction terms arise 
and cause the bounds on the norm to grow unboundedly as p ^ 00. 
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Theorem 4.1 (Stochastic Moser). For v > consider an invariant measure fx^ cind an associated 
stationary solution (x)^ of (3.2), where the stochastic forcing is assumed to be sufficiently smooth, e.g. (3.4) 
holds. Then the following bound holds 

E||w^||i<=o < C<oo (4.1) 

where C = C{a) is independent on v. 

An immediate consequence of estimate (4.1) is that any limit point fiQ of any sequence of invariant 
measures {^v}v>o is concentrated on i°°(T^). 

Theorem 4.2 (Invariant measure supported on L°°). Under the assumptions of Theorem 4.1, con- 
sider any collection of invariant measures {;Uv}v>o of (3.2). Then there exists a subsequence and a measure 
/io such that jiy. jiQ (weakly) in Pr{H^) as j — ' oo and jUo(L°°) = 1. 

Remark 4.3. Let us note that is an invariant measure for the Euler equation over L°° n X, where X 
is the fiber at r = of J^f. See Section 3 for details. 

We shall first give the proof of Theorem 4.2, assuming Theorem 4.1 holds, and then return and prove 
Theorem 4.1. 

Proof of Theorem 4.2. As in [Kuk04] by using (3.12) we have that {jUy} is tight, and hence weakly 
compact on Pr(i?°). Taking /Uq to be a limit point of {/iv}v>o in the weak topology of H^, there exists a 
sequence Vj such that 



for each 0eCfa(H°). 

According to Theorem 4. 1 we have 

sup E II w:^ II LOO = sup / WcJoWi'^diUvioio) < C < oo. (4.2) 

v>0 v>0 J 

We claim that this implies 

|wollL-rfA'o(i^o)<C. (4.3) 



Indeed, take to be a standard family of smooth moUifiers on U^. For i? > 0, £ > define (pR^^ito) = 
Wpe * a R. Young's inequality implies the (pR^e e ChiH^) so that 



^ (/)i^,£■(too)rf|Uo(wo) ^ C. 

Now, by Fatou's Lemma, we have 

(liminf llpc * cjoWl'^ a R)dLio[a)o] < C 

£->0 

Since ||toolli°° < liminff>o llpf * woIIl°° for each e L°° then (4.3) follows, completing the proof. □ 
Proof of Theorem 4.1. As a first step we rescale the time in (3.2). Taking T= t/v we obtain 

da) + 



^u-yS)-Aa)]dt = adW, 5(0) = ws(0), (4.4) 



V ) 

where we have denoted oj{t,x) = x), u{t,x) = u[t,x), and W(_t) = \/vWit). Note that Wit) has the same 
statistical properties as W[t). For ease of notation for we drop the tildes until (4.41) below. 

Fix T > 0, p > 1, and define Tjc to be an increasing sequence of times with To = and T^ as oo. Let 
Ik=lTk,2T]he a. sequence of time intervals approaching [T,2T]. 
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To analyze (4.4) we apply the Ito Lemma. This yields 



d\\CD\\1p 



(4.5) 



where we have denoted 



Ti„{t) = -\ u{t,x)-y(x){t,x)ci){t,x)\cD{t,x)\P ^dx = 
T2„{t)= I Aa){t,x)cD{t,x)\a)[t,x)\''~^dx 
73,p(t) = Ef Mt,x]\P-^a,„ixfdx 

Sm,pit)= I amix)(Oit,X)\(Dit,X)\P~^dx 

Jt2 



(4.6) 
(4.7) 
(4.8) 



for all p>2. In the identity for Ti p we have integrated by parts in x and used that V • m = 0. Since v > 0, 
we are dealing with spatially smooth solutions of (4.4), and hence the identity (4.5) may be justified by 
applying the Ito lemma pointwise in x, integrating over the torus, and using the stochastic Fubini theorem 
(see, e.g. [DPZ92]). Note however that (4.5) may be justified for much less spatially regular stochastic 
evolution equations, as recently established in [KrylO]. 

Let 5 e [Tfc, T^+i] and t> s. We start with (4.5) for p>2, integrated from s to t: 



mmlp -pf T2,p{j)dT = \\a){s)\\l, + E^R-Jl 
Js 2 

p pip-l) 

= I|w(s)||2p + 



[ T3,p{T)dT+pY,[ Sm,p{T)dW^' 
Js m Js 

[ T3,p{T)dT+pY,[ Sm,p{j)dW^' 
Js m JTk 



pL r S,n,p{T)dW,"'' 
m JTi- 



(4.9) 



where T2,p, T^^p, Sm,p are as defined in (4.6)-(4.8). We take the supremum of (4.9) over every t e Ik+i, and 
obtain 



0) 



p p{p-l) 



[ \T3,p{r)\dT + p sup W[ Sm.pMdW^" 
Jlk Kh+i I m Jit- 

+ p\L f Sm.pMdW,"" 

\ m JTu 

2 Jh t^Ii- \ m JTi- 



(4.10) 



where we have used that se [Tk, Tjc+i), and that T2,p < (cf. (4.18) below). This allowed us to bound from 
below the time integration on the left side from ls,2T] with the smaller one on [Tic+\,2T]. 

For the forthcoming computations it will be convenient to introduce the following standard notations for 
the stochastic (martingale) terms. For any < r < let 



Mir,t],p = Y.[' Sm.pMdW,"" 
m Jr 

and denote the running absolute maximum by 

M*:= sup W Sm.piT)dW," 

se[r,t] 1 m Jr 



(4.11) 



(4.12) 
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Finally the we define {M[r,-\p)t to be the quadratic variation of M[r,t],p and recall that (see e.g. [KS91]) 

{M^rAp)t= fY.Sm.pds. (4.13) 
Jr m 

Recall that by a version of the Burkholder-Davis-Gundy inequality given in [DMS05] (see also [RY99]) we 
have that, for any non-negative random variable Z, any r <t and any 5 > 

E(M* p V Zf < CBDGiS) ■ E«M[,,]p>J'2 V zf. (4.14) 

Here the constant Cbdg(5) is universal; it depends only on 5 and is independent of the form of the Martingale 
Mj* p or Z. Also, note carefully that there exists a 5o > such that 

Cbdg(5)<2^ whenever 5 < 5o. (4.15) 

This observation will be crucial below in estimates (4.33) and (4.37). 

We return to (4.10) and take an average of (4.10) for 5 e [Tfc, Tjc+i], and obtain 

ll^lll-(WP)-/' L T2,piT]dT 

1 rTk+i pfp — l^ r 

^^—f^ j^^ \\<ni,.Ln + ^^^ j^jr3,p(T)|rfT + 2pM[V,,2ri.p 

< + j^jr3,p(T)MT + 2pM(V„2^],p. (4.16) 

As usual in Moser iteration arguments, the lower bound on - T2,p is obtained by introducing v = \(o\f'^, 

2 

so that llwll^p = II f ll^2> and |Vf |^ = ^|Va)|^|a»|P~^. Then, upon integrating by parts in T2,p we have, point- 
wise in time, that 

-pT2,p = pip-l) [ \ya)\^\a)\P-^dx = 4^\\yv\\%>2\\yv\\% (4.17) 

for all p>2. Moreover, since we are in a two dimensional periodic box, the Sobolev embedding gives 

^lkl|2,. <||Vf;||2, + ||t;||2, (4.18) 

where 2* e [2,oo) is arbitrary, and the constant Cs > depends only on the the size of the box and the choice 
of 2* . Note that v is not zero mean in space and hence we need to add here a lower order term in (4. 1 8). Let 
us choose 2* = 4 for simplicity. Then, in view of (4.17) and (4.18), the left hand side of (4.16) is bounded 
from below as 

ll^lli~(/,,,;i^)+2||VHli.^,^,,.)>lklli.(,,^,,.)-2||.|li.(,^^^^,.)+2(|k||^ (4.19) 
By assuming that 

4|/fc+il=4(2r-rfc+i)<l, (4.20) 
which is automatically satisfied for all fc > if we ensure that 



r<i (4.21) 



we conclude from (4.19) that 



1 , 1 9 

>-\\v\\r^,r r2, + \\V\\,2,T r4v (4.22) 
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Let us now recall the following ij interpolation inequality. Suppose we have 1 < pi,p2, (^i,ci2, r\,r2< 
oo and < 7 < 1 satisfy 

1 r l-y 1 r 1— y 

- = — + -, - = — + '-. (4.23) 

r\ Pi qi r2 P2 qi 

Then, for any g e L^i (/; iP^) n I'?' (/; L*) we have 

llgllL'-i(/;L'-2) < \\S\\[piii-LP2)\\S\\\'h{i.Ln2) (4.24) 

with / c IR being some interval. Taking ri = 5, = 5/2, pi = oo, p2 = 2, = 2, = 4, and 7 = 3/5 in this 
inequality we find 

^ ||,.||2 ^ 11, .116/5 II ,,4/5 <i||i;||2 +— r4 25^ 

2^2/5 II -2c2/5ll^lli~(4+i;i2)ll^lll2(/hi;i*) " 2 " lli°=(4+i;i^) ^ 

by making use of the f- Young inequality. In summary, we have shown that the left hand side of (4.16) is 
bounded from below by 

^ c^"'^"l5p'2(/,+i;i=p'4) = ^ll'^lli2Ap(j^^j.iAp) (4-26) 

as long as (4.20) holds and for any p>2. Here have denoted = 2C^^ v 1, and denoted 

A=^. (4.27) 
For the term Ta ^-term on the left side of (4.16) we simply use Holder and obtain (pointwise in time) 

p{p-\) pip-V) p-2 2 

73,p< llwll^p ||o-||^p. 

Integrating the above on and using the Holder inequality in time we obtain 



r p(p-i) , p±? „-2 



(4.28) 



for all A; > 0. Thus from (4.16), (4.26), and (4.28), we obtain 



+ \\a\\l^+2pM*j.^^2np (4-29) 

with M*j^ p defined in (4.12), and we have used that WaWip ^ |T2|^'P||cr||i<x,. 
Let us now define 

K(p,r):=4C^ — -^ + \j^\2'Pp2ii^i2p +i+2p\J^\-p\h\'p\ (4.30) 

After some direct manipulations starting from (4.29), taking p^^ roots of both sides and then expectations 
we find that 

< Kip, T)pE[[\\aj\\L2Pu,;LP) V llo-|li»)'' V IT^fp |/fcr^M*j-^ 2r].p) ' 

<Kip,T)P CBDGip~'w((\Mi:^P(I,-LP)^ Ml'^Y ^ (4.31) 

which holds for all p > 2. Note that for the second inequality we used (4. 14). 



find 
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We next estimate the quadratic variation term, {M[Tf.,.],p}2j in (4.31). Starting from (4.8) and (4.13) we 



(1) 



I II ll''""^ 

llP|l'^ll^2p(7^.£p) 



Note that the second bound above makes use of the integral Minkowski inequality. 

Let us now summarize the estimates obtained, by combining (4.3 1) with (4.32). We have 

E(NllL2Ap(j^^^.iAp) V ||o-||ioo) <K{p,T)pCBDGiP~^W[\\^h^P(h;LP) V Ml""] (4.33) 

for all /? > 2. To set up a recurrence relation, it is hence natural to set p = pkin (4.33), where we define 

for all k>0, where we recall that A = 5/4. Let us now introduce some notation 

Ak = E(||w||i2p^;^._ip,) V ||c7||l<x>) 
J- _i 

ak = K{pk,T)PkCBDGiPic )• 



(4.32) 



Then, (4.33) reads 



So that 



E sup ||a»(f,-)llL<»<Ax)^ n^fc 



IE(llw||i4([o,2ri;L2) V ||o-||l°o). 



(4.34) 
(4.35) 

(4.36) 
(4.37) 

(4.38) 



te[r,2r] 

In view of (4. 15), we have that 

r\ 1^ logK{pk,Ty 
llflfc<Cexp 2^ 

fc>0 \k>0 Pic 

where C is a v- and T-independent constant. We now set 

rfc = rci-A"*'). 

Then, Tj,+i- = TX-''a + A'^) > rA-^/2= Tp^K We recall the definition of K[pk,T) from (4.30), which 
in view of the above choices may be bounded as 

K{pk,T)<4C's T--2pl + \J^\Pkpl{2T)'p^ +l+2pk\J^\pn2Typ^ <C/?2 T'l + i (4.39) 

where we have also used that < 2T < 1/4 (cf. (4.21)), and C is a sufficiently large T-independent 
constant. Using that Y.k>o P^^ = 5/2, and Y.k>o P^^logPfc < oo we may further obtain that 



n flfc <C{T-^'^ + 1)^'2 < c(r-s/4 ^ 



(4.40) 



fc>0 



for some sufficiently large v- and T-independent constant C. 

In summary from (4.37), (4.40) and recalling that these estimates were carried out for the rescaled 
equation (4.4) above we have in conclusion 



E sup ||w(?)||i<»<C(r ^''' + l)E(||w||i4([o,2r];L2)V||c7|lLoo), 
telT,2T] 



(4.41) 
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for any T < 1/8 and where C is independent of T and v. Rescaling to the original variable = S(vf) then 
with (4.41) we infer 

-r/v 



( C^ '^ \ 
E sup ||w(0||ieo<C(r"^"^ + l) / vU\\0J{s)\\n)ds 

tElTlv,2Tlv] I WO 



+ \\(J\\l'' 



(4.42) 



tE[Tlv,2Tlv] 

for any r< 1/8. 

We can now obtain the desired conclusion by taking o) to be a stationary solution of (3.2) corre- 
sponding to /iv- Recalling (3.14) and taking, for example T = 1/8 we have that 

E||w^||l»<E sup ||Ws(f)||i- 

t£[l/(8v),l/(4v)] 

<C N vEi\\col{s)\\l2)ds\ +||(7||l<. 



CM ^Eexp((5||w^(5)||2,)rfs 



1/4 



+ \\(J\\l- 



< Ca (4.43) 



for a constant independent of v. This gives (4. 1) concluding the proof of Theorem 4. 1 . □ 

Remark 4.4 (Linear drift-diffusion). Note that the L°° bounds obtained in Theorem 4.1 can also be 
shown to hold for any drift-diffusion equation 

d(ii) + [a-Vo)- Aa))dt = adW, V-a = 0, 

with sufficiently regular drift a and stochastic forcing a. Indeed as in (4.5) one may write the evolution of 
the LP-norm of o) with the analogous drift-term vanishing since a is divergence free. The rest of the proof 
follows without any change and one obtains that 

E sup ||w(?)||ic»<C(l + r-5/4)E[||w||i4([0 2r]-i2)V||c7||i.«), 

KlT,2T] 

for any < T < 1/8 and most importantly C is independent a. Note also that this estimate corresponds to 
the usual parabolic regularization in the deterministic case: weak-solutions are instantaneously in L°°. 

Remark 4.5 (Fractional Navier-Stokes). Note that the Moser iteration technique used to prove The- 
orem 4.1 may be used to obtain drift-independent L°° bounds for stationary solutions of the fractional drift- 
diffusion equation 

dcD + [a-ycD + [-A]'^'^a))dt = adW, 

for any power y e (0,2), where as in Remark 4.4 the drift a is divergence-free and sufficiently smooth. To 
see this, we recall the LP lower bound on the fractional Laplacian given in [CC04] 



which holds for any p>2. Using the 2D Sobolev embedding H^'^ c /^4/(2-7)^ ^j^g j-gpg^t fj^g argument 
given above in (4. 17)-(4.27), and obtain estimate (4.26) with Ay = 1 + j/A. Since for any y e (0, 2) we have 
A" — ' oo as n — ' oo the Moser iteration scheme may be completed mutatis-mutandis. In particular, setting 
a = ^u, which is divergence-free, in view of (4.4) one may use the above argument to study inviscid limits 
of the stochastic fractionally-dissipative Navier-Stokes equation. 

5. Modulus of continuity for the deterministic stationary problem 

In Section 4 we have proven that the stationary solution of (4.4) obeys v-independent bounds in L°°, 
that is E||aig||ioo is uniformly bounded in v. The key ingredients used in this argument were 

• Two-dimensionality: this ensures that the nonlinear term, whose size blows up (in comparison to 
the viscosity) as v — ► 0, vanishes altogether in estimates for the vorticity. To put it differently, 
there is no vorticity stretching term. 
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• Stationarity: this enables us to measure the L°° norm of the solution whose initial data is o)^ at 
time Ty~v~^, and hence obtain bounds on itself. 

Once we wish to estimate in more regular spaces, for example W with 5 > 1, or with y > 0, the 
nonlinear term does not vanish anymore, and since it's relative size becomes prohibitively large as v — ► 0, 
we do not seem to be able to obtain v-independent bounds on in spaces that are better than L°° (averaged 
over the probability space). 

In this section we exhibit an drift-independent bound, in a better norm than LP°, of solutions to the 
stationary drift-diffusion equation 

Lv = -\v + b-yv = f (5.1) 

for xeJ^, where b = b(_x) is a divergence free-vector field, but on which we have no bounds. The force is 
assumed to be in L°°, with zero-mean, and we consider solutions v such that fj2 vdx = 0. We view equation 
(5.1) as a deterministic toy-model describing the stationary solutions of (4.4) - the analogy is given by 
letting b = v~^u, v = 0), and noting that a){0) equals in law, for all f > 0. 

Following the ideas in [SSSZ12] in the spirit of [Leb07] we show that v obeys a logarithmic modulus 
of continuity which does not depend on the size of the drift b. In particular f is a uniformly continuous 
function. 

Theorem 5.1 (Modulus of continuity for deterministic stationary equation). Let b be divergence 
free and smooth, and v be a zero mean weak solution of {5.1 ), that is, ve and satisfies (5.1) in the sense 
of distributions. Then v obeys a drift-independent logarithmic modulus of continuity 

sup \v{x)-v[y)\< (5.2) 

\x-y\<r v^logl/r 

for some C > that is independent of b and all r e (0, r* ], for some universal constant r* . In particular, v is 
uniformly continuous. 

Proof of Theorem 5.1. Taking the inner product of (5.1), using the Poincare inequality and the fact 
that V • = 0, we obtain 

^1^11^2 + ^l|Vt;||2,<||/||^.||t;||i. 

which imphes 

||Vi;||2,<Cil|/||2, (5.3) 

for some Q > that is independent of b. We as usual denote oscjf v = maxfc - minjf v to be the oscillation 
of V over the set K. Fix xq = 0, Tq = diam(T^)/4A 1/2, and let By = Brixo) for any r < tq. Upon integrating in 
polar coordinates, dropping the normal derivatives, and using that by the ID Sobolev embedding [dBr) c 
C"idBr) for a e (0, 1/2) (see [SSSZ12, Theorem 4.2]), we obtain 

Ci\\f\\l,> dp, (5.4) 

Jr P 

for any r e (0, r^]. 

If we were able to establish that v is monotone in the sense of Lebesgue, i.e. to show that oscggp v is a. 
monotone function of p, the proof of the lemma would directly follow from (5.4). Instead we prove that v 
is almost monotone in the sense of Lebesgue, that is, up to an error of size r^. Let h solve 

Lh = -f in Br, h = on dBr 

so that L(_v + h) = in Br, and hence by the maximum principle 



OSCb,{V+ h) = OSCgB,.[V+ h) = OSCgBr V- 



(5.5) 
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We claim that 

osc5,h<C||/||i<-(5,.)r2 (5.6) 

for some constant C > that is independent of b. To prove (5.6), we rescale the problem to the unit ball by 
letting 

x= ry,h{_y) = hiry),b{_y) = bix),fiy) = /(x). 

It follows that 

-\yh + rb-Vyh = -r^f in Bi, h = OondBi. (5.7) 

We obtain the desired estimate by Moser iteration. Multiplying (5.7) by h\h\f~^ and integrating over Bi we 
obtain that for any p>2 

fc^f \vmP'')fdx={p-l)[ |V^|2|^|P-2rfx<r2|5i|i'P|l7llL=o|l^llj:;' 

P J Bi J Bi 

by using that V • = 0, and that h vanishes on dBi. Using the Sobolev embedding a in 2D, we get 

\\h\\l,^<Cpm^'PrHfMh\\l;^ (5.8) 

for some C > which is independent of p. Without loss of generality we take this constant C sufficiently 
large so that 2C\Bi\^'P > 1 for any p>2. Let 

Pfc = 2*^ and = max{||h||iPt, r^ll7llL<»} • 

It follows from (5.8) that 

ak,,<iCp,\B,\"P'=f'"a, 

and thus 

Whh^iB,) < C(r2||7||i«(B^) + 11^11^2(5^)) < Cr2 II 7|li»(B^) (5.9) 

for some C > 0. The last inequahty above follows by setting p = 2 in (5.8), and using the Holder inequality. 
Upon rescaling back to x-variables, (5.9) implies (5.6). 

Combining (5.5) with (5.6) gives that for any r < p < ro we have 

oscaBpt;>oscBpy-C||/||L-p^>oscB^y-C2ll/||i<»p^ (5.10) 
for some C2 > that is independent of b. Inserting this bound in (5.4) yields 

Ci 11/11^2 > ['^-(oscB,.i;-C2ll/|liocp2)2rfp. (5.11) 

Jr P 

We distinguish two cases, based on whether oscb,. v is larger or smaller than 20211/111°°/". When oscb^ v > 
2C2ll/|lL-r > 2C2ll/|lL-p2, then (5.6) implies that' 

p 



which implies 



On the other hand 



4^ll/lll2 



OSCb,. V < 



2C2I 



0SCB,.Z^<2C2ll/|lL-r< 



^logl/r 



for any r < < 1/2. The above two estimates imply (5.2). One may repeat this argument with xq being 
any point in T^, not just the origin, by periodically extending v and / to one more periodic cell, thereby 
concluding the proof. □ 
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In contrast, the parabolic case is more delicate. If we consider the linear problem 

dtv + bix, t)-Vv- Av = f, 

even if b is divergence-free, one may construct solutions that are not continuous functions for all time, 
although they obey the L°° maximum principle. See e.g. [SVZ12] for an example with rough drift. For the 
nonlinear problem 

one may hope to prove that in some average sense, the functions on the attractor remain continuous as v — 0. 
At the moment we do not know how to prove this. 

6. The Damped and Driven Navier-Stokes Equations and Other Scaling 

In this section we consider the weakly damped and driven stochastic Navier-Stokes equations 

du+iYu + u-Vu + yn-vAu]dt = v"pdW = v"Y,Pkdw'', V-u = 0, (6.1) 

k 

where Y = F^ y = tA"^' = t(-A)"^'^ and t > 0, y e [0, 1). As above for (3. 1) it follows immediately from the 
Kryloff-Bogoliouboff procedure that there exists an invariant measure p," for each v > 0, a eU. 
We now prove that a = is the only scaling of v in (6.1) which gives a nontrivial inviscid limit. 

Theorem 6.1 (Inviscid limits in different scalings). For a eU consider a collection of invariant 
measures {fJ,y}y>o of (6.1). Depending on the choice of a we have one of the following three scenarios: 

(i) If (X> 0, then for any vj — 0, we have that p". do, i.e. weakly in Pr{H^), where 5q is the Dirac 
measure concentrated at 0. 

(ii) lfa<0, then for any v„^0 such that the p"^ po then 

I \\ufj2dpo{u)=oo. (6.2) 

(iii) lfa = 0, then there exists a sequence Vn^O and a stationary martingale solution pa of 

du + {Yu + u-Vu + Vn)dt = pdW, V-m = 0, (6.3) 

such that Py^ pq. By a stationary Martingale solution of (6.3) corresponding to po we mean 
that there exists a stochastic basis 5fi := Wt)t>Qy^< W^) cind a predicable process 

M e ( [0, r] ; H") n l2 ( [0, r] , (6.4) 

which satisfies (6.3) and is stationary, i.e. the law Piuit) e A], Ac SS{H^), is independent of t 
and identically equal to po. 

Proof. Let be stationary solutions of (6.1) corresponding to p", and define to^ = V"*- • u^. Applying 
the Ito lemma to (6.1) and using stationaiity, we obtain: 

E(v||Vm^||2, + ||yi'2„v||2^-) ^ ^\\P\\12. (6.5) 



Additionally, by making use of the vorticity formulation of (6. 1), 

doj + iYcD + u-ycD-vAc!))dt = v"adW, a = y^p, (6.6) 
we also obtain, again with the Ito lemma and stationarity 



E(v||Vto-||2, + II Y^'^i^'Wl,) = ^llfTll^,. (6.7) 



Proof of (i). We begin with the case a > 0. From (6.7) we have 



y2a 

E||M^||^,_,„<— ||C7||2,. (6.8) 
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Using Chebyshev and compact embedding we infer that jj." is tight in Pr(//°). 

Consider any weakly convergent subsequence /u" ^ fi^. By the Skhorohod embedding theorem we may 
find a new probability space {n,^,P) and a sequence of H° valued random variables such that is 
equal in law to yu" and 

m'' ^ m" a.s. in H°, 
with uP equal in law to . Now with (6.8) and Fatou's Lemma we infer 

E||M°||r2 <liminf EllM^IIra <liminf ||o-||?2 = 0. 

^ v>0 ^ v>0 2t 

Hence = a.s. and therefore fi^ =6o- This proves the first item. 

Proof of (ii). Now we consider the case a < 0. For every v > we obtain 



vEIIVm^II^, = vE||w||^2 <CvE 



1^2 II 11^2 

<Cv^E(v||Vaj^||22 + II F^'^^^^ll^^) < v^+^^CIIall^^, 

where we have used that Y = tA~^, interpolation and the above balance relation (6.7). Combining (6.5) with 
the above estimate we obtain 

'^\\p\\l.<E\\Y'"u^\\l.+v^r-'"C\\a\\l., 

for constant C > which is independent of v. This implies 

v^"[\\p\\lj2-v^C\\(T\\l2)<E\\Y^'^u''\\l2<CE\\u''\\l2. (6.9) 
By assumption the assumption that p" p.^ as in (i), the Skhorohod embedding theorem yields a new 



probability space and which converges almost surely to m°, with the same laws as the original sequence. 
We infer 

limsupv2"(||p||^2/2-v^C||or||^2) <ClimsupE||M'' 11^2 <CE||m° 11^2, 

where we have used the Fatou lemma in the last estimate. This proves (6.2). 

Proof of (iii). Lastly we treat the case a = 0. By applying the Ito lemma to (6.6) we have that 

d\\a)''\\l, + {2\\Ya)''\\l2 + 2v\\a)''\\]j,)dt= \\p\\l,dt+ {p,cj''}dW 

Using stationarity we immediately obtain that 

If'' is uniformly bounded in L^in;L^{[0, T],H^-^'^)) (6.10) 

and moreover that 

M^(0) is uniformly bounded in L^(Q; H^~"^'^)) which implies that {Hv}v>o is tight on H°. (6.1 1) 
Returning to (6.1) and again making use of the Ito formula, 

dllw" 11^2 + (2||Fm'' 11^2 +2v||M''||^i)d?= \\a\\l2dt + {a,u''}dW, 
we infer with (6.1 1) that 

is uniformly bounded in L^in;L"°{lO, T],H°)). (6.12) 

In order to obtain a suitable compactness required to pass to the limit we need some additional uniform 
estimates on fractional the time derivatives of u^. We will apply the Aubin-Lions type compact embedding 

L^ilO, T],H'^~^'^) n W^'^'^ilO, T]; H~^) aci L^{[Q, T], H°) (6.13) 

and the Arzela-Ascoli type compact embedding 

1V^'^'*([0, T];H-^] cc C([0, T]; //"*), (6.14) 
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(see [FG95]). Define a sequence of measures {ju^}v>o on the path space C{[0, T];H~^] associated to {u^}v>o 
by 

H^{A) = P{u'' eA), A£SS{C[0,T];H~^]. (6.15) 
Using the embeddings, (6.13), (6.14), and suitable estimates we will next show that 

{fi^}v>o is tight in l2([o, T];H°) n C([0, T];H-^). (6.16) 
For this propose we write (6. 1) in its integral form 

u^m- I [Yu + Piu-yu)-vAu)ds +aW[t]:= lDit) + Isit), (6.17) 
V Jo I 

where P is the Leray projection operator onto divergence-free vector fields. Observe that, 

l|/D(f)lCi/3,4([o,ri;//-3)^C [^(||u^(0)||^_3 + IIFu^llH-3 + llvAu^||^_3 + IIP(M''-VM^)||^_3)^^f 

*J 

<C sup a + Wu^Wtf, (6.18) 
f£[o,r] 

where C is independent of 1 > v > 0, but may depend on T. By making use of a suitable version of the 
Burkholder-Davis-Gundy (see e.g [FG95]) inequality we have 

IE||/s(?)ll^i/3,4([o,r];H-3) - '^llo'lli2, (6.19) 
where again, by the assumption a = 0, C is independent of v but depends on T. Similar estimates yield 

l|/D(f)llL/3,2fro7-l-W-3)^'^ sup (l+llM^II^^), E||/s(f)||L/3,2f,0r,.H-3l^C'||O-||^2. (6.20) 

f£[o,r] u , J, 

To establish the first part of the tightness bound, (6. 16) we consider the sets 

:= {u : II M||i2([o,r],Hi-r/2) < i?} n {u : || m|| w^''^'H10,t];H-'^) - ^} 

so that, according to (6. 13) is compact in L^([0, T],H°) for every R> I. Observe that, with an appropriate 
application of Chebyshev's inequality and in view of (6.10), (6.12), (6.17), (6.20) 

V 1 C "^11 V \\2. ^ II II ^ ll2 

^ ((5^) hHlO,T],H^-ri^)~^^^»^Dit)\\w''^'HlO,T];H--') + ^^\\Isit)\\^^^^^ 

c c 

~^'^(""^"i^([0.ri,Hi-^'2)"^ ""^"i°°([0,T],H») + II^IIl^] ^ (6.21) 

where C is independent of v > and R>0. For the second half of the tightness bound (6.16) we define 

Bjf := {u: llM|lwi'W([o,r];H-^) --^l' 
and observe with (6.18), (6.17) (6.19) that 

2 8 

^^{[Bjf) ) <-E||/£,(?)||^vl/3,4([0,r];H-3)+ -^E||/s(f)||^l,3,4(jg .j,j.^_3) 

c c 

^;^E(||M-||2„„,^^,,,„, + ||a||4,)<-, (6.22) 

where, once again, C is independent of v > and R>0. With (6.21), (6.22) we may now infer (6. 16). 

With (6.16) in hand we now invoke the Skorokhod theorem obtain a sequence of processes {u^,W^) 
defined on a new probability space, {Q.,J^,P) such that 

u^^u almost surely in C([0,r];//"^)nI^([0,r];H°), (6.23) 
almost surely in C([0, r];ilo), (6.24) 
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is equal in law to and each W"^ is a cylindrical Brownian motion relative to the filtration given by 
the completion of a{iW^{s),u^(_s)) : s < t). Note that the sequence maintains the same uniform bounds 
as in (6.10), (6.12) and it follows that 

m" — M weakly in L^{n;L^ao,T];H^-^'^)), (6.25) 

u"^ u weakly* in L°°(O;i°°([0, T];H°)). (6.26) 

Finally, as in [FG95] that 5 is a stationary process on H° with stationary distribution ^q. 

An argument from [Ben95] may now be employed to show that, for each v > 0, solves (6.1), but 
relative to the new stochastic basis = {D.,^,P,^f,W^). Now by using the convergences (6.23)-(6.26) 
we may pass to the limit v — ► in (6.17) (with u^, W appropriately replaced with m^, W^) and establish that 
(m, W) satisfies (6.3) along with the required regularity. As in [FG95] it follows from the stationarity of 
and (6.23) that u is stationary. 
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